J.  Non-Newtonian  Fluid  Mech.  150  (2008)  104-1 15 


Journal of 
Non-Newtonian 
Fluid 

Mechanics 


www.elsevier.com/locate/jnnfm 


Homogeneous  flow  field  effect  on  the  control  of  Maxwell  materials 

Hong  Zhoua>*,  Wei  Kanga,  Arthur  Krenera,  Hongyun  Wangb 

a  Department  of  Applied  Mathematics,  Naval  Postgraduate  School,  833  Dyer  Road,  Bldg.  232,  SP-250,  Monterey,  CA  93943-5216,  United  States 
b  Department  of  Applied  Mathematics  and  Statistics,  University  of  California,  Santa  Cruz,  CA  95060,  United  States 

Received  16  August  2007;  received  in  revised  form  6  October  2007;  accepted  24  October  2007 


Abstract 

The  controllability  of  viscoelastic  fields  is  a  fundamental  concept  that  defines  some  essential  capabilities  and  limitations  of  the  resulting 
materials.  In  this  paper,  we  study  the  controllability  of  different  homogeneous  flow  fields  of  viscoelastic  fluids  governed  by  the  upper  convected 
Maxwell  model.  The  approach  is  largely  based  on  the  nonlinear  geometric  control  theory.  Through  the  analysis  of  the  control  Lie  algebra,  we 
find  the  submanifolds  in  the  state  space  on  which  the  homogeneous  flow  fields  are  weakly  controllable.  Our  approach  can  be  generalized  to  more 
complicated  systems. 
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1.  Introduction 

The  controllability  of  viscoelastic  fluids  has  important  applications  in  the  design  of  desired  materials.  An  important  issue  of 
controllability  is  the  question  whether  it  is  possible  to  steer  a  system  to  a  desirable  state  with  a  given  set  of  control  inputs. 

In  [5]  the  controllability  of  flows  of  linear  viscoelastic  fluids  was  investigated  by  Renardy  using  multi-mode  Maxwell  models. 
The  state  of  the  system  is  characterized  by  the  velocities  and  the  viscoelastic  stresses;  the  control  input  is  in  the  form  of  the  body 
force.  It  was  found  that  the  system  is  uncontrollable,  unless  the  initial  conditions  for  the  stresses  satisfy  a  set  of  constraints.  For  the 
degenerate  case  of  creeping  flow  where  the  density  is  zero  in  the  equation  of  motion,  there  is  no  controllability  unless  the  control 
input  is  distributed  along  the  entire  interval  in  the  physical  domain.  In  the  presence  of  inertia,  crucial  difference  occurs  between 
the  cases  of  one  or  more  relaxation  modes:  For  a  single  relaxation  mode  (i.e.  a  Maxwell  fluid),  exact  controllability  holds  provided 
the  time  interval  satisfies  certain  inequality;  for  multiple  relaxation  modes,  exact  controllability  holds  under  modified  regularity 
assumptions.  Another  piece  of  work  by  Renardy  [6]  was  focused  on  the  homogeneous  shear  flow  of  viscoelastic  fluids  with  several 
different  constitutive  models.  For  those  equations,  the  state  of  the  system  consists  of  viscoelastic  stresses  whereas  the  shear  rate  is 
regarded  as  a  control  input.  For  the  upper  convected  Maxwell  (UCM)  model,  it  was  revealed  in  [6]  that  the  reachable  set,  i.e.  the 
states  in  stress  space  which  are  accessible  from  a  given  initial  condition,  is  specified  by  a  positive  definiteness  inequality  of  the  stress 
tensor.  Very  recently  the  result  was  extended  to  the  control  of  nonhomogeneous  shear  flow  of  an  upper  convected  Maxwell  fluid  [7] 
where  the  state  of  the  system  and  the  available  control  are  the  same  as  those  in  [5]. 

The  goal  of  this  paper  is  to  extend  the  work  in  [6]  to  a  variety  of  homogeneous  flow  fields  on  the  control  of  viscoelastic  fluids.  Our 
work  differs  from  [6]  in  several  aspects.  First,  we  include  in  this  study  the  effect  of  different  homogeneous  flow  fields  in  addition  to 
the  shear  flow  addressed  in  [6].  Secondly,  due  to  its  nonlinear  nature,  our  approach  is  largely  based  on  nonlinear  geometric  control 
theory,  which  is  different  from  the  analysis  and  tools  used  in  [6]  for  the  UCM  model  under  shear  flow.  The  third  difference  lies  in  the 
definition  of  controllability.  Rather  than  the  general  concept  of  controllability,  in  this  paper  we  adopt  a  local  version  of  the  concept, 
namely  weak  controllability.  This  definition  has  been  widely  used  in  nonlinear  control  theory  as  well  as  in  engineering  applications. 
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In  many  real  life  applications  of  complicated  nonlinear  systems,  engineers  prefer  to  reach  a  distance  target  by  scheduling  a  sequence 
of  local  movements  for  the  reasons  of  model  uncertainties,  system  perturbations,  and  sensor  noise.  From  a  theoretical  viewpoint, 
Lie  brackets  of  vector  fields  provide  an  efficient  tool  to  prove  weak  controllability.  We  note  that  Lie  brackets  have  been  exploited  in 
[6]  for  a  nonlinear  model  and  the  Renardy  model.  But  the  systems  addressed  in  this  paper  are  different.  Among  the  various  types 
of  flows  addressed  in  this  paper,  the  special  case  of  shear  flow  was  also  addressed  in  [6] ;  and  the  results  are  consistent  except  that 
we  address  the  concept  of  weak  controllability,  not  the  reachable  set  derived  in  [6].  However,  for  the  case  of  extensional  flow,  we 
derive  the  constraints  of  the  reachable  set  in  Section  4. 

Similar  to  [6],  we  mostly  limit  our  study  to  the  two-dimensional  upper  convected  Maxwell  fluid  except  for  Section  5,  in  which  we 
introduce  some  remarks  and  discussions  on  three-dimensional  flows.  However,  similar  analysis  developed  in  this  paper  is  applicable 
to  higher  dimensions  (even  though  more  advanced  mathematical  tools  are  needed)  and  other  constitutive  models. 

We  organize  our  paper  as  follows.  First  we  give  a  brief  discussion  on  the  controllability  of  control  systems.  Then  we  introduce 
the  upper  convected  Maxwell  fluid.  After  that,  we  analyze  the  effect  of  various  homogeneous  flow  fields  on  the  controllability  of 
Maxwell  materials.  More  precisely,  extensional  flow,  shear  flow  and  general  planar  linear  flow  are  considered  where  the  submanifold 
of  controllability  is  characterized  for  each  case.  Then  we  provide  a  brief  discussion  on  3D  models.  Finally,  we  conclude  the  paper 
with  a  summary  of  the  main  results. 

2.  The  controllability  of  control  systems 

Consider  a  general  nonlinear  control  system  that  is  affine  in  control  [4] : 

m 

x  =  fix) + (i) 
(=1 

where  x  is  the  state  variable,  w*  e  M,  i  =  1, . . . ,  m,  are  the  control  variables.  In  general,  the  state  variable  takes  value  in  a  manifold 
of  dimension  n ,  denoted  by  M.  An  important  concept  for  systems  defined  by  (1)  is  its  controllability ,  which  characterizes  the  ability 
to  maneuver  the  system  from  one  state  to  other  by  proper  choice  of  control.  There  exists  a  huge  literature  on  the  controllability  of 
control  systems.  In  this  paper,  we  adopt  the  definition  and  follow  the  geometric  approach  from  [2]. 

Let  xo  be  a  point  in  M,  the  manifold  of  state  variables.  A  point  x\  in  M  is  said  to  be  reachable  from  xo  if  there  exist  piecewise 
continuous  input  functions,  wz-  =  oti(t),  so  that  the  trajectory,  x{t ),  of  (1)  with  initial  state  xo  reaches  x\  in  finite  time,  i.e.  x\  =  x(T) 
for  some  T  >  0.  For  nonlinear  control  systems,  the  global  reachability  is  usually  very  difficult  to  prove.  Instead,  a  practical  solution 
is  to  study  weak  controllability. 

Definition  1.  A  control  system  is  said  to  be  weakly  controllable  in  an  open  subset  D  c  M  if,  for  any  xo  e  D,  there  exists  an  open 
neighborhood  Uo  of  xo  so  that  the  set  of  points  reachable  from  xo  along  trajectories  inside  Uo  contains  at  least  an  open  subset  of  M. 

If  a  system  is  weakly  controllable,  it  implies  that  the  locally  reachable  states  form  a  “solid”  region.  More  restrictive  than  weak 
controllability,  the  concept  of  controllability  requires  that  any  two  points  in  M  are  reachable  from  each  other.  For  linear,  time  invariant 
control  systems,  controllable  and  weakly  controllable  are  equivalent.  However,  for  nonlinear  control  systems  like  those  studied  in 
this  paper,  the  determination  of  controllability  requires  the  global  geometric  properties  of  the  vector  fields  in  control  systems.  The 
description  of  the  set  of  points  reachable  from  a  given  point  is  still  an  open  problem  for  most  control  systems.  On  the  other  hand, 
weak  controllability  can  be  determined  by  using  the  dimension  of  the  control  Lie  algebra  and  a  distribution  generated  by  the  vector 
fields  associated  with  the  control  system;  nevertheless,  weak  controllability  implies  some  important  properties  of  a  control  system. 

In  (1),  /(x)  and  gfx)  are  vector  fields  on  the  manifold  M.  Under  the  Lie  bracket  operation,  [f  g],  the  space  of  smooth  vector 
fields  on  M  forms  a  Lie  algebra.  This  Lie  algebra,  the  smallest  subalgebra  containing  the  vector  fields/,  g\,  . . . ,  gm,  is  called  the 
control  Lie  algebra,  denoted  by  C.  At  each  point  x  e  M,  the  vectors  in  C  span  a  vector  space.  It  is  denoted  by  A <?(x),  i.e.: 

A<?(x)  =  span{X(x)|2f  is  a  vector  field  inC}. 

This  mapping  from  M  to  the  tangent  bundle  of  M  is  called  a  distribution.  The  following  is  a  useful  sufficient  condition  on  the  weak 
controllability  of  a  nonlinear  control  system. 

Definition  2.  A  control  system  satisfies  the  controllability  rank  condition  (CRC)  on  an  open  set  D  c  M  if 

dim(Ac(x))  =  n  (2) 

for  all  xeD,  where  n  is  the  dimension  of  the  manifold  M. 

Theorem  1  (Isidori  [2]).  A  control  system  of  the  form  (1)  is  weakly  controllable  on  an  open  set  D  if  it  satisfies  the  controllability 
rank  condition  on  D. 
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3.  The  upper  convected  Maxwell  (UCM)  model 


The  upper  convected  Maxwell  (UCM)  model  was  proposed  by  J.C.  Maxwell  over  a  century  ago.  It  has  been  widely  used  in 
polymer  rheology  for  the  description  of  a  viscoelastic  fluid  under  large  deformations  [1,3].  The  UCM  model  gives  a  viscoelastic 
constitutive  equation  and  can  be  written  in  the  form: 

T  -  (Vv)T  -  T(Vv)t  +  XT  =  2/iD,  (3) 


where  T  is  the  stress  tensor,  v  the  velocity,  Vv  the  velocity  gradient  tensor,  X  the  relaxation  rate,  /x  the  elastic  modulus  and  D  is  the 
rate-of-deformation  tensor  (i.e.  the  symmetric  part  of  the  velocity  gradient). 

We  consider  2D  homogeneous  viscoelastic  fluids  and  denote  the  stress  tensor  by 


T  = 


Tn 

Tn 


Tn 

T22 


(4) 


where  T\\  is  the  first  normal  stress  difference,  T22  the  second  normal  stress  difference  and  T\2  is  the  shear  stress.  Suppose  the  control 
input  is  denoted  y  and  it  is  which  is  closely  related  to  the  velocity.  Then  the  general  dynamic  problem  of  (3)  becomes 

T  =  F(y(0,  T),  T(0)  =  T0,  T(Wi)  =  Ti,  (5) 


where  To  and  Ti  are  the  given  initial  and  final  states,  respectively.  The  state  of  the  system  (5)  is  characterized  by  viscoelastic  stress 
T  with  three  components  T\\,  T22  and  T\2. 


4.  Flow  field  effect  on  the  controllability  of  Maxwell  materials 

Our  primary  purpose  in  this  section  is  to  study  the  effect  of  flow  fields  on  the  controllability  of  Maxwell  materials  for  various 
types  of  flows.  A  weak  controllability  condition  is  presented  for  each  case. 


4.1.  Extensional  flow 


For  a  fluid  in  a  homogeneous  extensional  flow  with  rate  y(t ),  the  velocity  is 

K-*«^ 


V  = 


so  the  velocity  gradient  is 


Vv  = 


yit) 

2 

0 


0 

m 


2  j 

The  rate-of-strain  tensor  becomes 
1  T  v(t) 

d  =  -[vv+(vv)T]=  Y~4y 

In  component  form,  (3)  becomes 

T\  1  -  (y(t)  ~  a)7’h  = 

Note  that  T\2  can  be  solved  exactly: 


1  0 
0  -1 


T\2  +  ^T\2  —  0,  T22  +  (y(t)  +  X)T22  —  —fiy(t). 


(6) 

(7) 

(8) 

(9) 


Tn  =  7i2(0)  exp (-Xt).  (10) 

Since  the  behavior  of  Tn  is  unaffected  by  the  control,  the  UCM  model  (3)  is  not  weakly  controllable  and  hence  does  not  satisfy  the 
CRC  anywhere.  However,  the  state  space  has  a  stable  and  control  invariant  subspace  Tn  =  0.  All  trajectories  of  the  system  under 
any  control  input  asymptotically  approach  this  subspace.  Therefore,  the  ultimate  behavior  of  the  control  system  is  represented  by 
the  reduced  system  on  the  stable  subspace: 


Tn  -( y(t)-X)Tn  =  fiy(t), 


T22  +  (y(0  +  ^)T22  —  —fi'Yit). 


(ID 
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For  mathematical  convenience,  we  introduce 


Xl 

>11' 

X2 

T22 

Then  the  system  (11)  can  be  rewritten  as 

dx  _ 

~T~  =  f(x)  +  g(x)u, 
at 

where 

U  =  y(t),  f(x)  = 

The  Lie  bracket  is 

[/.?]  =  VI-/-  v/.|  = 

Then  we  have 


X/JL 

—X/i 


—Xx\ 

11 +x  1 

— Xx2 

,  g(x)  = 

-IX  -x  2 

(12) 


(13) 


(14) 


(15) 


det[J,  [/,  g]]  =  det 


li +  x  i 

— /X  —  X2 


X/1 


— X ,  /i 


—  X/l(x  2  —  Xi). 


(16) 


It  follows  that  the  subsystem  (11)  of  the  UCM  model  under  extensional  flow  satisfies  the  CRC  and  hence  is  weakly  controllable 
on  the  set  of  all  states  where  the  two  normal  stress  differences  are  not  equal,  i.e.  X2  ^  x\  (or  T\\  /  T22). 

For  extensional  flows,  it  is  possible  to  characterize  the  set  of  reachable  states  for  the  subsystem  defined  by  (1 1)  which  is  equivalent 
to 


7n  =  —{X  -  y{t))T\\  +  /iy(t),  T2 2  =  -(X  +  y(t))T22  - 


Our  approach  below  is  inspired  by  Renardy’s  work  [6]. 
Introducing  new  unknown  functions: 

x  =  Tn  +/z,  y=T2  2  +  /x, 

the  ODE  system  (17)  becomes 

i;  =  —(X  -  y(t))x  +  X/jl,  y  =  —(X  +  y(t))y  +  X/jl. 

If  we  rescale  all  the  variables  as  follows: 


x  = 


/x 


y  = 


y_ 

/x’ 


t  =  Xt, 


m  = 


m 

X  ’ 


(17) 


(IB) 


then  after  dropping  all  tildes  for  simplicity,  (18)  has  the  form 

x  =  -(1  -  P(t))x  +1,  y  =  -(1  +  m)y  +  1.  (19) 

We  shall  prove  the  following  theorem. 

Theorem  2.  Let  x(t)  and  y(t)  be  solutions  of  the  system  (19)  with  initial  conditions  x(0)  and  y(0).  If  x( 0)  >  0,  then  x(t)  >  0  for 
t  >  0;  similarly ,  ify( 0)  >  0,  then  y(t)  >  0  for  t  >  0. 

Proof.  We  argue  by  contradiction.  Suppose  x(t)  =  0  for  some  r  >  0.  Let 

t\  =  inf{r|x(r)  =  0,  r  >  0}. 

Since  we  assume  that  there  exists  r  >  0  such  that  x(t)  =  0,  it  follows  that  t\  is  well  defined  and  is  finite.  Because  x(t)  is  a  continuous 
function  of  t ,  we  have  x(t\)  =  0.  Using  the  assumption  x(0)  >  0  and  using  the  definition  of  t\ ,  we  get  x(t)  >  0  for  0  <  t  <  t\.  Now 
we  evaluate  the  derivative  of  x(t)  at  t\  from  the  first  equation  of  (19)  to  obtain: 


=  =  l  >0. 

t=t\ 


dx 
d  t 


108 


H.  Zhou  et  al.  /  J.  Non-Newtonian  Fluid  Mech.  150  (2008)  104-115 


Thus,  we  have  x(t\  —  s)  <  0  for  s  positive  and  small  enough.  This  contradicts  with  the  result  obtained  earlier  that  x(t)  >  0  for 
0  <  t  <  t\ .  Therefore,  we  must  have  x(t)  >  0  for  t  >  0.  Similar  arguments  lead  to  the  conclusion  that  y(t)  remains  positive  for  t  >  0 
if  its  initial  value  is  positive. 

It  is  worthwhile  to  point  out  that  the  zero  stresses  T\\  =  T22  =  0  correspond  to  x  =  y  =  1. 

From  now  on,  we  shall  confine  our  attention  to  the  case  where  x(0)  >  0  and  y(0)  >  0.  To  find  the  reachable  set,  we  start  with 
two  lemmas.  □ 


Lemma  1. 

d  (xy) 
d  t 


The  solutions  x(t)  and  y(t)  of  the  system  (19)  satisfy 
=  x  +  y  —  2  xy. 


Proof.  The  proof  is  straightforward.  Using  (19),  we  have 
d(vy)  dx  dy 

—r—  =  —y  +  x-£  =  -(1  P(t))xy  +  y  -  (1  +  P(t))xy  +  x  =  x  +  y  -  2xy.  □ 
dt  dt  dt 


(20) 


Lemma  2.  The  solutions  x(t)  and  y(t)  of  the  system  (19)  satisfy 


d  In (y/x) 
dt 


1  1 

-2  P(t)+ - . 

y  x 


Proof.  Applying  the  chain  rule  and  (19),  we  find 

d  lny 


dlnx  ldr  1 

dt  x  dt  x 


dt 


Id  y  1 

= --77  =  -[!  +  £«] +  -. 

y  dt  y 


It  follows  immediately  that 

d  In  {y/x)  d(lny  —  Inx) 


1  1 

—  —2  P{t)  + - . 

y  x 


dt  dt 

Let  us  introduce  two  new  unknown  functions: 

y 


§  =  ^  =  ln\l~ 

Using  Eqs.  (20)  and  (21),  we  find 
d^xy  1  dixy) 


dt 

d  In  y/y/x 


1 

2^/xy  dt 
1  d  In  (y/x) 


1  1 

=(v  +  y  -  2xy)  =  - 
2Jxy  2 


X  + 


1  (i  i 

=  -P(t)  +  -[ - 

2  V37  x 


dt  2  dt 

In  terms  of  §  and  77,  the  above  ODE  system  becomes 

d§  dr)  1 

—  =  cosh  r)  —  §,  —  =  —  fit) - sinh  77. 

dt  d  t  i=  ' 


=  ~P(t) 


1  1 


Jxy2 


(21) 


(22) 


(23) 


Note  that  the  system  (23)  is  nonlinear.  Furthermore,  the  evolution  of  §  is  completely  determined  by  (§,  rj)  whereas  the  evolution  of 
T)  depends  on  the  control  parameter  Pit).  Of  course,  §  is  controlled  by  Pit)  via  its  effect  on  77.  □ 

Theorem  3.  For  t  >  0,  the  solution  %it)  of  (23)  satisfies 

m  -  1  >  e-'G(0)  -  1] 

with  equality  possible  only  if  rj(s)  =  0  throughout  the  interval  0  <  s  <  t. 

Proof.  From  the  first  equation  of  (23),  we  obtain 

m  - 1) 


dt 

which  gives 

d[e^ 


+  (§—!)  =  cosh  77  —  1 


1)] 


dt 


=  e^(cosh  77  —  1). 
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Integration  from  0  to  t  yields 

£(f)  -  1  =  e~*[£(0)  —  1]  +  e“r  [  e*[cosh  rj(s)  -  1]  ds. 

Jo 

Since  cosh  qis)  >  1,  it  follows  that 

$(0-  1  >e~'[£(0)-  1]  for?  >  0. 

The  equality  holds  only  when  cosh  qis)  =  1  for  0  <  s  <  t,  which  is  equivalent  to  qis )  =  0  for  0  <  s  <  t.  This  completes  the  proof 
of  Theorem  3. 

Let  (§o,  ho)  =  (§(0),  qiO))  denote  the  initial  condition  and  let  (§/,  ??/)  denote  the  state  we  would  like  to  reach  at  time  tf.  The 
reachable  set  at  time  t p  will  be  described  by  the  following  theorem.  □ 

Theorem  4. 

(1)  If  —  1  <  e-^  (§o  —  1),  then  (§/,  rjy )  is  not  reachable  at  time  tf. 

(2)  If%  f  —  1  =  e-^  (§o  —  1),  (§  /,  0  f)  is  reachable  at  time  tf  if  and  only  ifqo  =  qf  =  0. 

(3)  Ifl-f  —  l  >  (§o  —  1),  then  (§y,  r]  f)  is  reachable  at  time  tf. 

Proof. 


(1)  The  proof  follows  directly  from  Theorem  2. 

(2)  If  (£/,  r]f)  is  reachable  at  time  tf  and  —  1  =  (§o  —  1),  then  Theorem  2  implies  that  ^(s)  =  0  for  0  <  s  <  t.  Consequently, 

we  have  qi 0)  =  qitf)  =  0.  That  is,  qo  =  r)  f  =  0. 

On  the  other  hand,  if  q$  =  r\f  —  0,  we  select  fit)  =  0  to  keep  qit)  =  0.  This  leads  to  %(t)  —  1  =  e-?(§o  —  1)  and  thereby 
l;(tf)  =  §/.  Thus,  (§/,  q  f)  is  reachable  at  time  tf. 

(3)  The  proof  of  (3)  needs  a  lemma.  □ 


Lemma  3.  Consider  solving  the  ODE  system  (23)  forward  in  time.  Let  us  select  fit)  to  change  q  linearly  in  time  from  r\ 
t  =  0  to  q  =  q2  at  t  =  At.  Specifically ,  we  enforce  dri/dt  =  iq2  —  q\) / At  by  selecting 

1  Q2  —  Q] 

m  =  -W)  smh,(0-4f. 

For  this  choice  of  fit),  it  is  true  that 


q  i  at 


e-A/[£(0)  _i]  +  (i_  e_Af)(cosh^min  -  1)  <  £(A  t)  -  1  <  e-A?[|(0)  -  1]  +  (cosh^max  -  1)A t. 


where 


qmin  =  inf{|g||gis  between  q\  and ^2},  gmax  =  sup{|g||gis  between  q\  and^}- 


For  the  proof  of  this  lemma,  we  integrate  the  first  equation  of  (23)  to  obtain 

rAt 


f(A0-  1  =e-Af[§(0)-  l]+e“Ar 


/ 


e5[cosh  qis)  —  1]  ds. 


Since  qm{n  <  \rj(t)\  <  qm ax  for  0  <  t  <  At,  we  get  coshgmin  <  cosh  qis)  <  coshgmax.  Thus,  it  follows  that 

e“Ar[?(0)  -  1]  +  (1  -  e-Af)(cosh^min  -  1)  <  $(Af)  -  1  <  e~At[m  ~  1]  +  (cosh^max  -  1)A t. 

Now  we  turn  our  attention  to  the  proof  of  part  (3).  We  first  select  fit)  to  change  q  linearly  in  time  from  q  =  qo  at  t  =  0  to  q  =  q  at 
t  =  At.  Then  we  choose  fit)  to  keep  q  at  q  =  q  from  t  =  At  to  t  =  tf  —  At.  Finally,  we  select  fit)  to  change  q  linearly  in  time  from 
q  =  q  at  t  =  tf  —  At  to  q  =  qyatt  =  tf.  However,  the  condition  t-(tf)  =  l=f  will  not  be  automatically  satisfied.  Let  us  consider  how 
to  choose  At  and  q  to  make  i=(tf)  =  £/. 

For  the  special  case  of  q  =  0,  we  look  at  the  upper  bound  on  §(£/).  Applying  Lemma  3  for  each  of  the  three  time  sub-intervals 
[0,  At],  [At,  tf  —  At]  and  [tf  —  At,  tf],  we  have 

f(A t)  -  1  <  e-Af(fo  -  1)  +  (cosh^o  -  l)Af,  §(f/  —  At)  —  1  =  e-(f/-2Af>[£(A t)  -  1], 

%(tf)  -  1  <  e_Af[£0/  -  At)  -  1]  +  (cosh  rif  -  1)A t. 
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Combining  these  results,  we  obtain 

%(tf)  -  1  <  -  At)  -  1]  +  (cosh T)f  -  1)A t  <  [£(Af)  -  1]  +  (cosh??/  -  l)Ar  <  e_f/(§0  -  1) 

+  (cosh  r\ o  —  l)At  +  (cosh  Of  —  1)A t. 


Recall  that  in  part  (3)  we  assume  §/  —  1  >  e~tf  (§o  —  1).  So  we  can  select  At  small  enough  such  that 

%(tf)  —  1  <  %f  —  1  forg  =  0.  (24) 

Once  such  a  At  is  found,  we  fix  it.  For  the  general  case,  we  look  at  the  lower  bound  on  §(£/).  Applying  Lemma  3  for  each  of  the 
three  time  sub-intervals  [0,  At],  [At,  tf  —  At]  and  [tf  —  At,  tf],  we  get 

$(A 0  -  1  >  e-Af(£0  -  1),  Rtf  -  At)  -  1  >  e-(f/-2A°K(A0  -  1]  +  [1  -  e-(f/-2Af)](coshg  -  1), 

$(*/)  -  1  >  e~At[i;(tf  -  At)  -  1], 

Putting  all  these  results  together  yields 

%(tf)  -  1  >  z~At[$(tf  -  At)  -  1]  >  er(tf~At\^{At)  -  1]  +  e_Ar[l  -  e“(?/“2A,)](cosh^  -  1) 

>  e“?/(§ o  -  1)  +  e“Ar[l  -  e“(r/“2Ar)](coshtf  -  1). 


Let  us  select  q  large  enough  such  that 

%{tf)  —  1  >  §/  —  1  for  ^sufficiently  large.  (25) 

Since  §(£/)  is  continuous  in  q ,  it  follows  from  (24)  and  (25)  that  we  can  find  a  value  of  q  such  that  %(tf)  =  £/.  This  concludes  the 
proof  of  Theorem  4. 

In  summary,  the  reachable  set  is  R  +  Ro  where  R  and  Ro  in  the  (§,  ??)  plane  are 

7?  =  {(I/,  ??/)||/  -  1  >  e-^o  -  1)},  Ro  =  {($/,  »?/)l$/  -  1  =  e_f/(?0  -  1)  and??/  =  W  =  0). 


In  the  (.v,  y)  plane,  /?  and  /?o  are  given  by 

R  =  {(*/,  yf)\sNJy~f  -  1  >  e_f/(v'x0y0  -  1)},  *o  =  {(*/,  yf)W*Jyf  -  1  =  e~tf(Vxm  -  1)  and  ^  =  —  =  1}. 

X  f  X0 

In  terms  of  (T\\ ,  T22)  and  the  original  time  before  scaling,  R  and  Ro  have  the  expressions: 

1  > 


R=  <(Tn(tf),T22(tf)) 


Tn(tf) 


+ 


OP?-) 


7n(0) 


+  1 


722(0) 


+  1-1 


7?o 


(7’n(r/),722(t/)) 


722(0) 


+  1-1 


andTn(O)  =  r22(0),  rn(r/)  =  r22(f/) 


4.2.  57z£<2r /tow 

For  a  shear  flow  with  rate  y(0,  the  velocity  is 

v  =  (y(0y,  0).  (26) 

The  rate-of- strain  tensor  becomes 

1  T  y(t) 

D=  -[Vv  +  (Vv)T]  =  ^ 

The  UCM  model  (3)  becomes 

Til  -  2y(t)T\2  +  XTn  =  0,  T\2  -  y(t)T22  +  k T\2  =  tiyit),  722  +  kT22  =  0.  (28) 


0  1 
1  0 


(27) 


H.  Zhou  et  al.  /  J.  Non-Newtonian  Fluid  Mech.  150  (2008)  104-115 


111 


Then  722(0  =  722(0)  exp(— Xt)  so  the  system  is  not  weakly  controllable.  We  consider  the  subsystem  where  T22  =  0: 

T\\  -  2y(t)T\2  +  XTn  =  0,  T\2  +  XTn  =  liyit). 

The  system  (29)  can  be  expressed  as 

dx  _ 

—  =  f(x)  +  g(x)u, 

where 


Xl 

>11" 

*2 

Tn  _ 

The  Lie  bracket  is 
0 


[/.  = 


So 


X/jL 


u  =  y(t). 


fix) 


— Xx\ 
-XX2 


g(x) 


2X2 

IX 


det[g,  [/,  g]]  =  2Xhx2. 


(29) 

(30) 

(31) 

(32) 

(33) 


Therefore,  the  subsystem  (29)  of  the  UCM  model  under  shear  flow  satisfies  the  CRC  when  the  shear  stress  is  nonzero,  i.e.  X2  ^  0 
(or  T]2  /  0)  and  hence  is  weakly  controllable. 

By  manipulating  the  system  (29),  it  was  found  that  the  set  of  reachable  states  of  (29)  is  given  precisely  by  the  following  inequality 

[6]: 

fiTu(tf)  -  Tl2(tf )2  >  er^fUiTniO)  -  T12(0)2].  (34) 

Note  that  the  condition  (34)  is  stronger  than  our  result  in  this  special  case  because  it  gives  global  reachable  sets  whereas  ours  is  a 
condition  on  local  controllability  only. 


4.3.  General  planar  linear  flow 


Now  we  consider  a  general  planar  linear  flow  with  velocity 
v  =  (p\(t)x  +  p2(t)y,  P3(t)x  -  p\{t)y). 


(35) 


This  general  case  is  different  from  the  previous  ones  because  an  invariant  subspace  may  not  exist.  As  a  result,  we  have  to  deal  with 
the  full  system  model  rather  than  a  reduced  control  system.  The  rate-of- strain  tensor  is 


D=-[Vv  +  (Vv)A]  = 


p\(t) 

P2(t )  +  P3(0 
2 


P2(t)  +  P3(t)  - 
2 


(36) 


Then  the  UCM  model  (3)  becomes 


T\\  -  2{p\{t)T\\  +  p2(t)T\2)  +  XTn  =  2p,pi(t),  T\2  -  (p2(t)T22  +  P2,{t)T\\)  +  XT\2  =  /lipiit)  +  Pi(t)), 

f22  -  2(p3(t)Tl2  -  pi(t)T22)  +  XT22  =  -2ppi(t). 

The  system  (37)  can  be  cast  into 

-7-  =  fix)  +  g\ix)u\  +  g2ix)u2  +  glix)U3, 

where 


X  = 

Xl 

x2 

— 

' Tn ' 
Tl2 

,  Ml  =  Plit),  U2  =  P2it),  M3  =  P3(t),  f(x)  = 

1 

1  1 
>* 

X  X 

to 

_ 1 

X3 

T22 

—XX3 

2pt  +  2x\ 

2X2 

0 

0 

,  hX)  = 

/X  +  X3 

-  hX)  = 

n  +  X] 

— 2/x  —  2x3 

0 

2x2 

(37) 


(38) 


(39) 


It  is  practical  to  consider  the  situations  where  only  one  pft)  is  a  control  parameter  and  the  other  two  are  constants.  Therefore,  we 
consider  three  cases. 
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Case  1.  p\{t)  is  the  only  control  parameter  and  other  parameters  are  constants. 
The  system  (37)  can  be  rewritten  as 

dx  -*• 

777  =  fl(x)  +  gi(x)ui, 
at 

where  x,  g\ ,  u\  are  defined  in  (39)  and 


fl(x)  = 


—Xxi  +  2X2P2 

—Xx2  +  (/X  +  X^)P2  +  (/x  +  X\)P2 
-XX2  +  2x2  P3 


After  some  calculations  we  obtain 


[hdfuhUfulfuh]]] 


2/i  T  2x\ 
0 

— 2/x  —  2^3 


4x2p2  +  2Xp  %P2Pi{pt  +  Xi)  +  2X2p, 

-2p3(lA  +  Xl)  +  2p2(tl  +  X3)  4Xfl( P2  -  P3) 

-4X2P3-2X/1  -&P2P3(lX  +  X3)  -  2X2p 


(40) 


(41) 


(42) 


and 

det([gi,  [fuhl  (fi,  t/i,|i]]])  =  8A2/i(xi  -  X3)(x\p3  -X3P2)  +  SXfi(p2  -  Pi) 

x  [Xp(x i  -  x3)  -  4p,x2(p2  ~  Pi)  +  4x2(xip3  -  X3p2)\  .  (43) 


•  If  p2  =  p3,  then  the  determinant  (43)  is  not  zero  if  x\  /  x3,  which  implies  that  the  CRC  is  satisfied  and  the  UCM  system  (37)  is 
weakly  controllable  there. 

•  If  p2  /  P3,  we  set  the  determinant  (43)  to  zero  and  solve  for  x2  in  terms  of  x\  and  x3.  After  some  algebras,  we  find  that 

Kx\  ~  X3)[p-(p2  -  Pi)  +  {X\P3  -  X3p2)] 

x2  = -  (44) 

4(p2  -  Pi)[p(pi  -  Pi)  -  (xiPi  -  X3P2)) 

provided  that  n(p2  —  pi)  —  (xi  p$  —  x3 p2 )  /  0.  In  other  words,  if  p2  /  p3  and  //( p2  —  Pi)  —  (x\ p2  —  x2p2)  /  0,  then  the 
determinant  (43)  does  not  vanish  if  x2  does  not  lie  on  the  surface  described  by  (44).  As  a  result,  the  system  (37)  satisfies  the  CRC 
and  is  weakly  controllable  there. 

•  If  p2  7^  Pi  and  pi p2  —  pi)  —  (x  1  pi  —  Xip2)  =  0,  then  the  determinant  (43)  simplifies  to 

detail,  [/1,  ia,  [/1,  t/i,  it]]])  =  i6kVcp2  -  wx*  i  -  *3).  (45) 

So  the  system  (37)  satisfies  the  CRC  and  is  weakly  controllable  if  x\  ^  X3. 


To  summarize,  we  consider  the  state-parameter  space: 
R5  =  {(xi,  X2,  Xi,  P2,  Pi)\Xi,  Pj  €  R}. 

Define  some  surfaces  in  R 5  by 


■S)  =  {{x\,X2,X3,  p2,  Pi)\p2  =  Pi),  S2  =  {{xi,X2,Xi,p2,Pi)\p{p2- Pi)-(x\Pi  —  Xip2)  —  0}, 


s3 


|(Xl,X2,X3,  P2,  Pi) \x2 


XQci  -  Xj)[p,(p2  -  Pi)  +  (MPi  -  XjpX]  \ 

4(p2  -  Pi)[^(P2  -  Pi)  -  (x\Pi  -  XiP2.)]  )  ’ 


s4  =  {(X\,X2,Xi,  p2,  Pi)\x\  =Xi). 


(46) 


Then  in  the  state-parameter  space,  R5,  the  system  is  weakly  controllable  in  R5  \  (S)  U  S2  U  S3}.  In  Si,  the  system  is  weakly 
controllable  in  5)  \  At;  In  S2,  the  system  is  weakly  controllable  in  S2  \  {Si  U  S4}. 

Case  2.  p2(t)  is  the  only  control  parameter. 

The  system  (37)  can  be  rewritten  as 

dx  -*•  ^  ^ 

77  =  fl(x)  +  g2(x)u2, 
at 


(47) 
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where  x ,  g2,  U2  are  defined  in  (39)  and 


fl(x)  = 


-Xx\  +  2{p  +  x\)p\ 
-1X2  +  O  +X\)pj, 

-A.X3  -  20  +  xh)P\  +  2X2P3 


Then  we  have 


[g2,  Ui,  gi\ ,  \h-  [  fi ,  g2]]]  = 


and 


2x2  2(xi  +  p)p3  —  4x2/2i  -4p\  P3(n  +  X] )  +  2kjip3  +  8x2/22 

/i+X3  —2p\{p  +  X3)  +  k/i  4/2jO  +  X3)  -  2.p\(p  +  *1)  +  A./i(A.  -  2p\) 

0  -2/23(/i  +  X3)  4piP3(p  +  X3)-4pi(X2P3  +  Xp) 


det([g2.  [?2,  I2],  [/2,  l.?2, 12]]])  =  4A/x/?3  /?302  -  *3)  +  2/23x1  (/X  +  x3)  +  A.x2(x3  -  p)  -  2/23X2 j  . 
If  we  set  the  determinant  (50)  to  be  zero  and  solve  xi  in  terms  of  x2  and  X3,  we  find  that 


_  P3 (2X2  +  *3  _  M2)  _  XX2(x3  -  jl) 

Xl  2/230  +  X3) 

provided  that  the  denominator  /t  +  X3  /  0.  From  this,  we  have  the  following  conclusions. 


(48) 


(49) 


(50) 


(51) 


•  If  p  +  X3  /  0  and  x  does  not  lie  on  the  surface  given  by  (5 1 ),  then  the  determinant  (50)  is  not  zero,  the  CRC  holds  and  the  system 
(37)  is  weakly  controllable. 

•  If  p  +  X3  =  0,  i.e.  X3  =  —p,  then  the  determinant  is  reduced  to 

det([g2,  [/2,  g2],  [/2,  [/2,  g2]]])  =  —8 Xpp3x2(Xp  +  P3X2).  (52) 

Hence  the  determinant  does  not  vanish  if  x2  /  0  and  x2  7  —Xp/p3.  In  this  case,  the  CRC  holds  and  the  system  (37)  is  weakly 
controllable. 


Let  us  summarize  our  results  geometrically.  In  R 5  =  {(xi,  x2,  X3,  p\,  P3)\xi,  pj  e  R],  dehne  three  surfaces  given  by 

{»3(2x?  +  X?  —  it2)  —  A.X2(X3  —  it)  I 
(x\,X2,X3,  P1,P3)\X\  =  - - - 7 - ; - 7 -  >  , 

2/73(/X  +  JC3)  J 

S7  =  {(X1,X2,X3,  pi,  P3)\X2  =  0  OYXfl  +  P3X2  =  0}.  (53) 

Then  the  system  is  weakly  controllable  at  all  points  in  R5  \  {S5  U  S^}.  In  S5,  the  system  is  weakly  controllable  at  all  points  in  S5  \  S7. 

Case  3.  p^it)  is  the  only  control  parameter. 

The  system  (37)  can  be  rewritten  as 

dx  -*•  _ 

-7-  =  Mx)  +  g3(x)U3,  (54) 

at 

where  x,  g3,  w3  are  defined  in  (39)  and 


%(x)  = 


-Ax  I  +  2(/i  +  x\)p\  +  2  X2P2 

-XX2  +  (p  +  X3)/22 

-XX3  -  2(/i  +  X3)p\ 


Then  we  have 


[?3.  [?3,  ibL  t/3,  [?3,  k3]]] 


0 

P  +  X] 

2x2 


— 2(xi  +  P)P2 

2p\{p  T  xi)  +  Xp 
2p2(p+X3)  +  4x2pi 


4/21  /22(/t  +  Xl)  -  4/22(x2/22  +  A/i) 

4/22(/r  +  xi)  -  2p\{p  +  X3)  +  Xp(X  +  2p\) 
4/2i  /22(/r  +  X3)  +  8/22x2  +  2p2Xp 


and 


det([|3,  l.?3, 13].  1%,  1/3, 13]]])  =  ~4Xpp2 


p2(p2  -  x2)  +  2 p2X3(p  +  xi)  +  kx2(xi  -  p)  -  2 p2x\ 


(55) 


(56) 


(57) 
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As  before,  we  set  the  determinant  (57)  to  be  zero  and  solve  X3  in  terms  of  x\  and  X2  and  obtain 

piili  -  x?)  +  Xx2(x\  -  n)-  2p2xl 
=  - - - —  ^ 

-2p2(hL+Xl) 

if  /x  +  x\  /  0.  Several  observations  immediately  follow: 


(58) 


•  fi  +  x  1^0  and  X3  does  not  lie  on  the  surface  (58),  then  the  determinant  (57)  is  nonzero,  the  CRC  holds  and  the  system  (37)  is 
weakly  controllable. 

•  If  /z  +  x\  =  0,  i.e.  x\  =  —pi,  then  the  determinant  (57)  becomes 


det(g3,  [/3,  £3],  t/3,  t/3,  £3]]])  =  8X[ip2X2(X[l  +  /?2*2), 

which  is  not  zero  when  X2  /  0  and  X2  7^  —Xpi/p2.  Then  the  CRC  holds  and  the  system  (37)  is  weakly  controllable. 
To  summarize,  consider  three  surfaces: 


(59) 


S&  =  {(X\,X2,X3,  p\,  P2)\pb+X\  =0},  S9=  <  (X\,X2,X3,  pi,  P2)\X3 


s  10  =  {(*l,*2,*3,  PI,  P2)\X2  =  0  or  +  p2X2  =  0} 


P2(p  ~  x\)  +  XX2(X\  ~  pi)  ~  2p2xl  1 
-2j92(/X+Vl)  J 


(60) 


in  the  state-parameter  space  R5  =  {(xi ,  JC2,  JC3,  p\ ,  J92)} •  Then  the  system  is  weakly  controllable  at  all  points  in  R5  \  {Ss  U  5*9}.  In 
5g,  the  system  is  weakly  controllable  at  all  points  in  S%  \  5io- 

5.  Some  remarks  on  3D  models 


It  should  be  pointed  out  that  we  study  the  2D  fluids  mainly  due  to  their  mathematical  convenience.  In  fact,  if  we  assume  that  the 
stress  tensor  is  of  the  following  form 

T11  T\2  0 

T\2  T22  0 
0  0  0 

then  the  3D  UCM  model  is  reduced  to  the  2D  model  considered  in  this  paper.  If  one  consider  the  full  3D  UCM  model  coupled 
with  a  three-dimensional  flow,  then  the  analysis  will  be  much  more  complicated.  As  an  example,  consider  the  3D  UCM  model  in 
the  presence  of  a  3D  elongational  flow  where  the  velocity  field  is  described  by  v  =  (—y(t)(x/2)  —  y(t)(y/ 2),  y(t)z)  and  the  UCM 
model  (3)  can  be  expressed  as 


dx  -> 

-7-  =  f(x)  +  g(x)u, 
at 


where 


x  = 


"xi" 

Til' 

X2 

Tn 

X3 

Tl3 

X4 

T22 

X5 

T23 

_  X6  _ 

_  T32  _ 

The  Lie  brackets  are 
'  —  Xpi 
0 
0 

Xpi 
0 

_  2Xpi  _ 


[/,£]  = 


=  g  u 


u  =  y(t),  f(x)  =  —Xx, 


[/,  5il  =  ^g  =  g  2,  •  •  ■ 


(61) 


g(x)  = 


M  -xi 
~X2 

1 

2'X3 

—  pi  ~  X4 

1 

2X5 

2  jx  -|-  2x6 


(62) 


(63) 
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It  is  obvious  that  the  determinant  of  the  Lie  brackets  vanishes  and  one  cannot  conclude  weak  controllability  in  3D  case.  However, 
it  is  always  possible  to  find  control  invariant  submanifolds  in  the  state  space  that  are  weakly  controllable.  The  characterization  of 
such  invariant  submanifold  needs  geometric  analysis  [2]  which  is  beyond  the  scope  of  this  paper. 

Generally  speaking,  for  3D  models  coupled  with  3D  flows,  the  controllability  becomes  more  subtle  and  complicated.  One  needs 
to  use  geometric  approach  to  identify  controllable  submanifolds. 

6.  Concluding  remarks 

In  this  paper  we  have  applied  the  controllability  rank  condition  to  the  vector  fields  in  the  upper  convected  Maxwell  model  to 
study  the  controllability  of  viscoelastic  fluids  driven  by  various  homogeneous  flow  fields.  In  our  control  system,  the  state  variable  is 
the  stress  and  the  available  control  is  related  to  different  flow  rate.  For  the  upper  convected  Maxwell  model  coupled  with  different 
flow  fields,  we  find  that 

•  The  UCM  model  under  extensional  flow  is  weakly  controllable  on  a  stable  invariant  submanifold  when  the  first  normal  stress 
difference  T\\  does  not  equal  to  the  second  normal  stress  difference  T22.  The  reachable  set  is  also  derived. 

•  The  UCM  model  under  shear  flow  is  weakly  controllable  on  a  stable  invariant  submanifold  when  the  shear  stress  T\2  is  nonzero. 

•  The  UCM  model  under  general  planar  linear  flow  is  weakly  controllable  in  some  subsets  of  the  state  space.  The  constraints  that 
define  these  subsets  are  characterized  explicitly  by  equations  of  the  stress  components. 
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